TECHNICAL  REPORT 


Defense  Threat  Reduction  Agency 
8725  John  J.  Kingman  Road,  MS-6201 
Fort  Belvoir,  VA  22060-6201 


DTRA-TR-15-16 


Development  of  Mouse  Lung 
Deposition  Models 

Distribution  Statement  A.  Approved  for  public  release;  distribution  is 
unlimited. 


July  2015 


DTRA01-03-D-0014 


Dr.  Bahman  Asgharian 
Prepared  by: 

Applied  Research  Associates 
8537  Six  Forks  Road 
Suite  600 

Raleigh,  NC  27615 


DESTRUCTION  NOTICE 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  to  sender. 


PLEASE  NOTIFY  THE  DEEENSE  THREAT  REDUCTION 
AGENCY,  ATTN:  DTRIAC/  J9STT,  8725  JOHN  J.  KINGMAN  ROAD, 
MS-6201,  ET  BELVOIR,  VA  22060-6201,  IE  YOUR  ADDRESS 
IS  INCORRECT,  IE  YOU  WISH  THAT  IT  BE  DELETED  EROM  THE 
DISTRIBUTION  LIST,  OR  IE  THE  ADDRESSEE  IS  NO 
LONGER  EMPLOYED  BY  YOUR  ORGANIZATION. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

00-07-2015  Technical  February  10,  201 1  -  October  7,  2012 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 


Development  of  Mouse  Lung  Deposition  Models 


DTRA01-03-D-001 4-0030 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Dr.  Bahman  Asgharian 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Applied  Research  Associates,  Inc. 

8537  Six  Forks  Rd. 

Suite  600 

Raleigh,  NC  27615 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Defense  Threat  Reduction  Agency 
8725  John  J.  Kingman  Road,  MS  6201 
Fort  Belvoir,  VA  22060-6201 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


Distribution  Statement  A.  Approved  for  public  release;  distribution  is  unlimited. 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

DTRA 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

DTRA-TR-15-16 


14.  ABSTRACT 

Mathematical  models  were  developed  to  predict  deposition  of  inhaled  particles  in  the  lungs  of  two  strains  of  mice.  The 
deposition  model  included  an  algorithm  to  calculate  particle  inhalability.  Model  predictions  included  filtering  capabilities  in  the 
nasal  passages,  and  deposition  in  various  regions  of  the  lung.  Particle  inhalability  was  smallest  in  mice  compared  with  that  for 
larger  rodents  and  humans.  In  addition,  there  was  a  significant  deposition  of  particles  in  the  nasal  passages,  with  the 
consequence  of  limited  penetration  of  particles  into  the  lung.  As  a  result,  particle  deposition  in  lung  airways  was  small, 
particularly  for  larger  particles.  The  findings  from  the  study  indicated  that  mice  inhalation  studies  should  be  conducted  with 
sub-micrometer  particles  to  allow  for  sufficient  deposition  in  the  lung  to  produce  the  desired  biological  responses. 


15.  SUBJECT  TERMS 


Mouse  Lung  Geometry  Lung  Deposition  Modeling  Inhalability  Regional  Deposition 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

Unclassified 


b.  ABSTRACT 

Unclassified 


c.  THIS  PAGE 

Unclassified 


17.  LIMITATION 

18.  NUMBER 

OF  ABSTRACT 

OF  PAGES 

BAR 

28 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Richard  Fry 


19b.  TELEPHONE  NUMBER  (include  area 
code) 

703-767-3193 


standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


CONVERSION  TABLE 

Conversion  Factors  for  U.S.  Customary  to  metric  (SI)  units  of  measurement. 


MULT  I  PL  Y- 
TO  GET  ◄ - 


->  BY 


->  TO  GET 


BY  ◄- 


DIVIDE 


angstrom 

1.000  000  X  E  -10 

meters  (m) 

atmosphere  (normal) 

1.013  25  X  E  +2 

kilo  pascal  (kPa) 

bar 

1.000  000  X  E  +2 

kilo  pascal  (kPa) 

barn 

1.000  000  X  E  -28 

meter^  (m^) 

British  thermal  unit  (thermochemical) 

1.054  350  X  E  +3 

joule  (J) 

calorie  (thermochemical) 

4.184  000 

joule  (J) 

cal  (thermochemical/cm^) 

4.184  000  X  E  -2 

mega  j  oule/m^  (MJ/m^) 

curie 

3.700  000  X  E  +1 

*giga  bacquerel  (GBq) 

degree  (angle) 

1.745  329  X  E  -2 

radian  (rad) 

degree  Fahrenheit 

tfc  =  (t°f  +  459.67)71.8 

degree  kelvin  (K) 

electron  volt 

1.602  19  X  E  -19 

joule  (J) 

erg 

1.000  000  X  E  -7 

joule  (J) 

erg/ second 

1.000  000  X  E  -7 

watt  (W) 

foot 

3.048  000  X  E  -1 

meter  (m) 

foot-pound-force 

1.355  818 

joule  (J) 

gallon  (U.S.  liquid) 

3.785  412  X  E  -3 

meter^  (m^) 

inch 

2.540  000  X  E  -2 

meter  (m) 

jerk 

1.000  000  X  E  +9 

joule  (J) 

joule/kilogram  (J/kg)  radiation  dose 

absorbed 

1.000  000 

Gray  (Gy) 

kilotons 

4.183 

tera joules 

kip  (1000  Ibf) 

4.448  222  x  E  +3 

newton  (N) 

kip/inch^  (ksi) 

6.894  757  x  E  +3 

kilo  pascal  (kPa) 

ktap 

1.000  000  X  E  +2 

newton-second/m^  (N-s/m^) 

micron 

1.000  000  X  E  -6 

meter  (m) 

mil 

2.540  000  X  E  -5 

meter  (m) 

mile  (international) 

1.609  344  X  E  +3 

meter  (m) 

ounce 

2.834  952  x  E  -2 

kilogram  (kg) 

pound-force  (lbs  avoirdupois) 

4.448  222 

newton  (N) 

pound-force  inch 

1.129  848  X  E  -1 

newton-meter  (N-m) 

pound-force/inch 

1.751  268  X  E  +2 

newton/meter  (N/m) 

pound-force/foot^ 

4.788  026  X  E  -2 

kilo  pascal  (kPa) 

pound-force/inch^  (psi) 

6.894  757 

kilo  pascal  (kPa) 

pound-mass  (Ibm  avoirdupois) 

4.535  924  x  E  -1 

kilogram  (kg) 

pound-mass-foot^  (moment  of  inertia) 

4.214  011  X  E  -2 

kilogram-meter^  (kg-m^) 

pound-mas  s / f oot^ 

1.601  846  X  E  +1 

kilogram-meter^  (kg/m^) 

rad  (radiation  dose  absorbed) 

1.000  000  X  E  -2 

**Gray  (Gy) 

roentgen 

2.579  760  X  E  -4 

coulomb/kilogram  (C/kg) 

shake 

1.000  000  X  E  -8 

second  ( s ) 

slug 

1.459  390  X  E  +1 

kilogram  (kg) 

torr  (mm  Hg,  0°  C) 

1.333  22  X  E  -1 

kilo  pascal  (kPa) 

*The  bacquerel  (Bq)  is  the  SI  unit  of  radioactivity;  I  Bq  =  I  event/s. 
**The  Gray  (GY)  is  the  SI  unit  of  absorbed  radiation. 


TABLE  OF  CONTENTS 


CONVERSION  TABLE . iii 

LIST  OE  EIGURES . v 

PREEACE . vi 

SUMMARY . vii 

I.O  INTRODUCTION . I 

1. 1  MOTIVATION . I 

2.0  APPROACH . 3 

2.1  INHALABILITY  MODELING . 3 

2.2  MODELS  OE  URT  EILTERING  EEEICIENCY  BY  IMPACTION  AND  BROWNIAN 

DIEEUSION . 5 

2.3  LUNG  GEOMETRY . 8 

2.4  LUNG  VENTILATION  AND  PHYSIOLOGY  PARAMETERS . 8 

2.5  MECHANISTIC  MODELOE  PARTICLE  DEPOSITION  IN  THE  LUNG . 9 

2.6  SOETWARE  IMPLEMENTATION . 9 

3 .0  MODEL  PREDICTIONS . II 

4.0  CONCLUSIONS . 15 

5.0  REEERENCES . 17 

6.0  DEEINITIONS,  ACRONYMS,  AND  ABBREVIATIONS . 19 


LIST  OF  FIGURES 


Figure  2-1.  Compartmental  representation  of  the  lung  geometry . 6 

Figure  3-1 .  Partiele  deposition  in  the  lungs  of  miee  via  endotraeheal  breathing . 1 1 

Figure  3-2.  Deposition  fraetion  of  various  size  partieles  at  different  lung  depths . 12 

Figure  3-3.  Partiele  deposition  in  the  lungs  of  B6C3F1  and  BALB/e  miee  via  nasal  breathing.  13 


PREFACE 


The  research  work  described  in  this  report  was  conducted  for  the  Defense  Threat  Reduction 
Agency  (DTRA)  under  contract  number  DTRA01-03-D-0014-0030.  The  Contract  Officer’s 
Representative  for  this  effort  was  Mr.  Rick  Fry  of  DTRA’ s  Information  Systems  Capability 
Development  Office  (J9CBI). 


SUMMARY 


Assessment  of  the  deposited  dose  in  the  lung  enables  linking  of  exposure  eoneentration  with 
biological  responses.  Various  strains  of  mice  are  used  in  inhalation  studies  to  study  toxicological 
end  points.  However,  little  or  no  information  is  available  on  particle  deposition  in  the  lungs  of 
mice.  Lung  deposition  models  were  developed  for  particle  inhalation  in  the  respiratory  tract  of 
BALB/c  and  B6C3F1  mice.  The  deposition  model  included  particle  inhalability  and  nasal  losses. 
Particle  inhalability  in  mice  was  lower  than  that  in  rats.  In  contrast,  deposition  of  the  same  size 
particle  was  higher  in  mice  nasal  passages  than  that  in  rats.  Thus,  fewer  particles  entered  the 
mouse  lung  in  comparison  with  rat  particle  inhalation.  The  penetration  was  severely  limited  for 
micrometer  and  larger  particle  sizes.  Therefore,  relevant  and  not  identical  particle  sizes  to  human 
exposure  scenarios  must  be  used  for  interspecies  data  extrapolation.  Size  relevancy  issues  are 
addressed  by  considering  compatible  deposition  of  particles  in  mice  and  humans. 


1.0  INTRODUCTION 


The  use  of  various  strains  of  miee  in  toxieology  studies  has  been  on  the  rise  in  reeent  years  (e.g., 
[1-5])  due  to  the  relatively  low  eost  of  breeding  and  husbandry  eompared  with  other  speeies, 
short  life  eyele  to  observe  long-term  health  effeets,  and  appropriateness  of  miee  as  the  surrogate 
model  for  eertain  human  biologieal  end  points.  Under  a  federally  funded  researeh  projeet  by  the 
National  Institutes  of  Health,  deposition  models  were  developed  for  two  strains  of  miee  (B6C3F1 
and  BALB/e).  These  models  were  ineluded  in  our  multiple-path,  partiele  dosimetry  model 
(MPPD  V2.1;  Applied  Researeh  Assoeiates,  Ine.,  Raleigh,  NC)  to  prediet  the  deposited  dose  of 
inhaled  biologieal,  ehemieal,  and  radiologieal  agents  under  various  exposure  seenarios.  Coupling 
of  mouse  deposition  models  with  those  of  humans  and  other  speeies  in  MPPD  allows  for  dose 
extrapolation  as  well  as  injury  and  easualty  risk  analysis  for  seenarios  under  whieh 
measurements  are  available  in  one  speeies  but  not  the  others.  The  deposition  model  was  based  on 
detailed  meehanies  of  airflow  and  partiele  transport  in  the  nasal  passages  and  lung  airways.  The 
small  geometrie  dimensions  of  the  nostril  and  low  minute  ventilation  limited  the  penetration  of 
airborne  partieles  into  the  respiratory  traet  of  miee.  Thus,  models  of  partiele  inhalability  for 
miee  were  also  developed  and  integrated  into  MPPD.  Model  details  and  implementation  in 
MPPD  are  deseribed  below. 

1.1  MOTIVATION 

Injury  and  easualties  resulting  from  inhalation  of  an  agent  are  direetly  related  to  the  deposited 
dose  in  the  nasal  or  lung  airways.  In  miee,  inter-strain  variability  in  the  respiratory  traet 
geometry,  breathing  parameters,  and  innate  sensitivity  warrant  strain-speeifie  investigation  into 
deposition  and  ensuing  biologieal  outeome.  The  availability  of  measurements  of  the  deposited 
dose  in  the  lungs  of  miee  is  limited.  In  addition,  experiments  are  expensive  and  time-eonsuming, 
and  eannot  eover  all  exposure  seenarios.  Alternatively,  as  for  other  speeies,  meehanistie  models 
ean  be  developed  to  assess  internal  dose  in  miee  for  various  exposure  and  breathing  seenarios. 
Dosimetry  models  provide  an  inexpensive,  quiek,  and  effieient  means  to  evaluate  the  injury  or 
easualty  risk  assoeiated  with  inhaled  exposure  to  an  agent. 

Meehanistieally-based  partiele  deposition  models  were  developed  to  deseribe  the  transport  and 
deposition  of  inhaled  partieles  in  the  lungs  of  two  strains  of  miee  (B6C3F1  and  BALB/e).  The 
deposition  models  were  based  on  the  multiple-path  geometry  that  was  previously  developed  for 
humans,  rats,  and  rhesus  monkeys  [6],  [7].  Inputs  to  the  model  ineluded  lung  geometry  and 
volumes,  and  physiologieal  parameters  of  miee.  In  addition,  models  for  inhalability  and  nasal 
deposition  of  partieles  were  developed  based  on  information  for  rats  and  ineluded  in  MPPD  for 
miee. 
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2.0  APPROACH 


The  geometry  of  the  respiratory  tract  is  fairly  complex.  It  is  made  up  of  thousands  of  multi¬ 
scaled  airways  with  poorly  defined  shape.  In  addition,  there  is  a  lack  of  information  on  airway 
dimensions  for  the  entire  respiratory  tract,  which  includes  the  tracheobronchial  and  pulmonary 
regions.  Furthermore,  airflow  and  particle  transport  through  the  respiratory  system  is  not  quite 
fully  understood.  Thus,  detailed  modeling  of  particle  transport  in  the  entire  respiratory  system  is 
not  possible  with  the  given  limited  information.  However,  there  have  been  simplifications  for 
various  aspects  of  the  respiratory  tract  geometry  and  airflow  and  particle  transport  based  on 
observations  and  experimentations.  These  simplifications  have  been  incorporated  into  the 
deposition  model.  Consequently,  development  of  a  respiratory  tract  deposition  model  is  divided 
into  several  steps,  which  describe  the  holistic  process  of  particle  transport  through  lung  airways 
from  which  the  deposited  dose  of  inhaled  particles  can  be  computed.  The  steps  describe  the 
derivations  of  models  for  particle  inhalability  and  upper  respiratory  tract  (URT)  deposition  of 
particles,  creation  of  mice-specific  lung  geometries  and  ventilation,  and  finally,  mathematical 
formulations  for  particle  penetration  and  deposition  through  lung  airways.  While  some  steps  are 
species-invariant  such  as  ventilation,  others  are  not  (e.g.,  inhalability  and  URT  deposition 
calculations).  Details  of  each  step  are  presented  in  the  following  sections. 

2.1  INHALABILITY  MODELING 

As  the  particle  size  becomes  larger,  its  increasing  inertia  begins  to  hamper  its  ability  to  follow 
the  inhaled  airflow  streamlines.  If  the  particle  deviates  significantly  from  the  airflow  streamlines, 
it  may  fail  to  enter  the  URT  through  the  nostrils.  Considering  a  suspension  of  airborne  particles 
being  inhaled  into  the  respiratory  tract,  the  ability  of  particles  to  enter  the  URT,  or  their 
inhalability,  decreases  with  increasing  particle  size.  Therefore,  models  for  particle  inhalability 
must  be  developed  to  correct  for  the  reduction  of  exposure  concentration  due  to  escaping 
particles.  Particle  inhalability  is  lower  in  mice  than  larger  size  rodents  because  of  the  smaller 
nasal  openings  of  mice.  Particles  of  inhalable  sizes  in  rats  and  larger  rodents  may  not  necessarily 
be  inhalable  by  mice.  Thus,  the  characterization  of  particle  inhalability  is  more  significant  for 
mice  than  other  rodents. 

Semi-empirical  relationships  have  been  constructed  for  particle  inhalability  as  a  function  of 
particle  inertia  based  on  either  laboratory  measurements  of  particle  penetration  through  the  URT 
or  computational  fluid  dynamics  analysis,  which  solves  the  airflow  field  and  particle  transport  in 
the  neighborhood  of  the  URT  during  inhalation.  However,  no  such  studies  have  so  far  been 
conducted  in  mice  and  a  major  data  gap  exists  which  prevents  direct  development  of  particle 
inhalability  models  in  mice.  Considering  that  mice  and  rats  have  similar  nasal  structures, 
differing  primarily  in  size  and  breathing  rates,  one  can  perform  dimensional  analysis  to  derive 
equivalent  inhalability  expressions  in  mice  from  existing  models  in  rats.  While  particle 
inhalability  is  related  to  particle  inertia,  a  literature  review  reveals  that,  with  the  exception  of  one 
study  [8],  particle  inhalability  has  exclusively  been  related  to  its  diameter.  Hence,  the  effect  of 
breathing  rates  (and  hence  inertia)  was  not  accounted  for  in  these  studies.  Inhalability  models 
that  depend  on  particle  diameter  alone  provide  first-order  approximation  and  are  not  used  here. 
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We  adopted  the  inhalability  model  of  Asgharian  et  al.  [8]  for  Long-Evans  rats  and  extended  it  to 
miee  based  on  dynamie  similarities  between  rats  and  miee. 

Asgharian  et  al.  [8]  eonducted  a  series  of  short-term  nose-only  inhalation  exposure  studies  in 
whieh  Long-Evans  rats  were  exposed  to  monodisperse  aerosols  ranging  in  diameter  from  0.9  pm 
to  4.2  pm.  Undisturbed  aerosol  eoneentration  measurements  in  the  exposure  environment  and 
measurement-based  ealeulations  of  URT  and  LRT  deposition  fraetions  were  subsequently  used 
to  eonstruet  a  semi-empirieal  relationship  for  partiele  inhalability,  or  inhalable  fraction  (IE),  in 
rats  as  a  function  of  the  product  of  particle  aerodynamic  diameter  ( d^g )  squared  and  inhalation 

flow  rate  ( Q ). 


IF  =  1- 


(1) 


where  d^gQ  is  the  Stokes  parameter  and  is  given  in  g  x  pm^  /  s  with  d^g  and  Q  denoting  the 

aerodynamic  diameter  and  inhalation  flow  rate,  respectively.  Coefficients  a=19.87  and  P= 
-0.7466  were  found  by  fitting  equation  (1)  to  the  calculated  inhalability  values  [8].  It  is  assumed 
that  equation  (I)  equally  applies  to  mice  as  well,  because  rats  and  mice  have  similar  facial  and 
URT  geometries,  that  is  inhalability  equations  in  rats  and  mice  are  related  to  the  impaction 
parameter  (Stokes  number,  Stk)  by  the  same  relationship.  Thus, 


Stk  = 


dLQ 

18pDh 


Rats 


18pDj 


Mice 


(2) 


in  which  p  is  the  absolute  viscosity  of  air  and  D,^  is  the  hydraulic  diameter  of  the  URT. 
Hydraulic  diameter  is  defined  with  respect  to  the  nasal  openings  in  each  species  because  the 
Stokes  number  indicates  the  inertial  strength  of  particles  to  enter  or  bypass  the  nasal  opening.  By 
simplifying  and  rearranging  equation  (2),  the  following  results  are  obtained. 


D 


h-Rat 


D 


xdaeQ 


h-Mouse 


Mouse  ' 


(3) 


which  can  be  substituted  in  equation  (1)  to  find 


IF  =  1 
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1-1- a 


D 


h-Rat 


V  ^h-Mouse  J 


Mouse 


(4) 


where  =  1.065  mm  and  =  0.42  mm  are  hydraulic  diameters  for  the  rat  and 

mouse  nasal  openings,  respectively.  Equation  (4)  is  further  simplified  by  substituting  for 
hydraulic  diameters. 
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IF  =  1 


1 


(5) 


Where  a'=  9.92.  Inspection  of  equation  (5)  reveals  that  inhalability  in  mice  is  lower  than  that  in 
rats  for  the  same  Stokes  parameter.  When  inhalability  predictions  in  mice  and  rats  are  compared, 
there  is  an  additional  reduction  of  inhalability  in  mice  due  to  reduced  minute  ventilation  and, 
hence,  a  lower  Stokes  parameter. 

2.2  MODELS  OF  URT  FILTERING  EFFICIENCY  BY  IMPACTION  AND 
BROWNIAN  DIFFUSION 

Deposition  of  particles  in  the  mouse  URT  occurs  by  two  distinct  mechanisms.  For  sub¬ 
micrometer  and  smaller  particles.  Brownian  diffusion  is  responsible  for  the  removal  of  particles 
from  the  inhaled  air.  In  contrast,  inertial  impaction  is  the  dominant  mechanism  of  deposition  of 
particles  in  the  nasal  passages  for  sub -micrometer  and  larger  particle  sizes.  Particle  deposition  by 
sedimentation  is  negligible  in  the  mouse  URT  due  to  the  small  residence  time  of  particles  in  the 
airway  passages. 

Approaches  for  calculating  particle  deposition  in  the  URT  of  mice  are  similar  to  those  proposed 
for  inhalability  assessment,  which  are  either  based  on  experimental  measurements  or 
computational  fluid  dynamics  analysis.  In  this  case,  neither  approach  is  feasible  due  to  the  lack 
of  information.  Therefore,  the  models  developed  in  other  rodents  with  similar  URT  geometry  to 
mice  will  be  extended  to  mice  by  making  proper  adjustments. 

There  are  no  computational  fluid  dynamics  studies  available  for  deposition  of  sub-micrometer 
and  larger  particles  sizes  in  the  URT  of  mice.  However,  particle  deposition  has  been  measured  in 
the  URT  and  LRT  of  CDl  mice  [9].  Raabe  et  al.  [9]  conducted  a  nose-only  exposure,  in  which 
CDl  mice  were  exposed  briefly  to  monodisperse  aerosols  ranging  in  size  from  0.27  pm  to  9.65 
pm.  The  exposure  was  composed  of  mainly  fine  and  coarse  particles.  Hence,  the  reported 
measurements  reflected  deposition  primarily  by  inertial  impaction  and  included  deposited  masses 
of  particles  in  the  nasal  and  lung  airways.  Exposure  atmosphere  was  also  sampled  at  an  available 
exposure  port  of  the  nose-only  exposure  unit.  These  measurements  were  used  here  to  calculate 
regional  deposition  fraction  of  particles  by  dividing  the  deposited  mass  in  a  given  region  by  the 
mass  of  inhaled  particles  during  the  inhalation  study.  In  a  simple-minded  scenario,  one  may 
assume  the  respiratory  tract  is  represented  by  URT  and  LRT  compartments  while  breathing  is 
approximated  to  be  steady-state,  justifiable  by  the  fact  that  the  Womersley  number  (ratio  of 
inertial  to  viscous  forces  -  a  measure  of  flow  unsteadiness)  is  very  small.  It  may  further  be 
assumed  with  reasonable  accuracy  that  URT  deposition  efficiencies  during  inhalation  and 
exhalation  are  similar.  Referring  to  Figure  2-1,  deposition  fractions  in  the  head  and  lung  are 
found  from  the  following  relationships: 
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Figure  2-1.  Compartmental  representation  of  the  lung  geometry, 

UluRT  =  IHin  ^  'HuRT-imp  +  (^^in  “  ^  'HuRT-imp  “  ^^LRt)  ^  'HuRT-imp  (6) 

in  which  riuRT-imp  the  deposition  effieiency  of  particles  in  the  URT  by  impaction  during 
inhalation  and  exhalation,  is  the  mass  of  inhaled  partieles,  andm^Rj  and  mLRj  are  masses  of 
particles  deposited  in  the  URT  and  LRT,  respectively.  Dividing  both  sides  of  equation  (6)  by 
mj„  and  solving  for  gives 
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In  whieh  Aup,.j.  =  muj^.j.  /m^^  and  A^rx  =  m^Rj  /  mj^^  are  deposition  fractions  of  particles  in  the 
URT  and  LRT.  Measured  mass  of  deposited  partieles  in  the  URT  and  LRT  from  the  experiment 
of  Raabe  et  al.  [9]  were  used  to  ealeulated  regional  deposition  fractions  and  URT  deposition 
efficieney  (equation  (7))  for  eaeh  partiele  size  used  in  the  experiment.  Caleulated  URT 
efficieneies  were  fit  to  a  logistie  function  suggested  by  Zhang  and  Yu  [10]  to  obtain  an 
expression  for  particle  deposition  by  impaetion  in  the  URT  of  mice. 
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Sinee  animals’  minute  ventilations  were  not  reported  by  Raabe  et  al.  [9],  the  empirieally-derived 
equation  of  Guyton  [11]  was  used  in  equation  (8)  to  find  a  =  3.96  and  P  =  0.1 17  with  R^  =  0.98 . 

There  is  no  information  on  deposition  of  ultrafine  partieles  in  the  URT  of  mice  either  by 
measurements  or  theoretieal  modeling.  Comparison  of  the  nasal  structure  of  mice  with  the  same 
for  rats  reveals  similar  dominant  features  while  exhibiting  typical  fine  intra-  and  inter-subjeet 
variations.  Thus,  for  predictions  with  order-of-magnitude  level  of  aeeuraey,  the  nasal  geometries 
of  the  two  speeies  are  assumed  to  be  similar,  which  allows  deposition  models  for  ultrafine 
partieles  in  rats  to  be  extended  to  miee.  Based  on  measurements  in  the  nasal  casts  of  rats,  Cheng 
et  al.  [12]  obtained  the  following  expression  for  losses  of  ultrafine  partieles  in  the  nasal  passages 
of  rats  by  Brownian  diffusion  during  inhalation  and  exhalation. 


fi  URT-diff 
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where  D  is  the  partiele  diffusion  eoeffieient  ( em^  /  see  )  and  Q  is  the  inhalation  or  exhalation 
flow  rate  (L/min).  Coeffieients  a,  P,  and  y  were  found  by  eurve  fitting  the  measurements  to 
equation  (9)  [12],  While  deposition  effieieney  by  Brownian  diffusion  in  equation  (9)  is  explieitly 
deseribed  in  terms  of  the  diffusion  eoeffieient  and  breathing  flow  rate,  the  more  general  form  of 
the  equation,  whieh  will  inelude  inter-subjeet  variability,  is  given  by 

77^  (10) 


where  Sehmidt  number  (Se)  and  Peelet  (Pe)  numbers  are  given  by 
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In  whieh  v  is  the  kinematie  viseosity  of  air.  Substituting  for  So  and  Pe  in  equation  (10)  and 
rearranging  variables  gives 


1-e 


(12) 


Beoause  of  geometrio  similarity  of  nasal  passages,  equation  (12)  denotes  deposition  effieieney  by 
diffusion  in  the  nasal  passages  of  both  rats  and  mioe.  Equating  the  exponential  term  of  equation 
(12)  for  the  ease  of  the  two  speoies  gives 
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By  letting  p  =  -(P'+y')  and  y  =  y' ,  it  beoomes  apparent  that  deposition  effieieney  by  diffusion  in 
the  URT  of  mioe  oan  be  deseribed  by 


hd  =l-e 


Dh_ 

Dh_ 


Rats 

Mice 


dP  O^ 

^Mice^Mice 


(14) 


Henoe  equation  (9)  oan  be  used  to  oaloulate  URT  losses  in  mioe  simply  by  sealing  the  eoeffieient 
a  by(D[j_j^j^^g  /Dh-mice)^  •  H  1^  noted  that  hydraulio  diameters  in  equation  (14)  are  defined  based 
on  volume  and  surfaoe  area  of  the  nasal  passages  whereas  those  in  equation  (4)  are  based  on  the 
perimeter  and  oross-seotional  area  of  the  nostril.  Replaoing  for  the  values  of  hydraulio  diameters 
in  rats  and  mioe  it  is  found  that  a  =  13.688  and  1 1.407  for  inhalation  and  exhalation, 
respeotively,  while  P  =  0.517  and  y  =  -0.234  remain  the  same  as  for  rats. 
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2,3  LUNG  GEOMETRY 


Measurements  of  lung  airway  parameters  sueh  as  length,  diameter,  gravity  and  bifureation  angles 
are  difficult  to  make  due  to  the  enormous  number  of  airways  and  complexity  of  airway  shape  and 
lung  structure.  Measurements  are  typically  made  along  selected  pathways  spanning  from  the 
trachea  down  to  terminal  bronchioles  for  the  conducting  airways.  The  same  measurements  are 
made  for  the  pulmonary  airways  including  respiratory  bronchioles,  alveolar  ducts,  and  alveolar 
sacs  in  the  selected  paths  and  are  complemented  by  measuring  the  size  of  the  alveoli  and  the 
number  of  alveoli  per  duct.  Modeling  each  airway  as  a  cylindrical  tube  and  the  lung  as  a 
branching  network  of  airways,  representative,  typical-path,  dichotomous  and  asymmetric  lung 
geometries  have  been  constructed  for  a  variety  of  species. 

14-generation,  typical-path,  symmetric  tracheobronchial  airway  geometries  for  two  strains  of 
mice  (B6C3F1  and  BALB/c)  were  constructed.  The  method  proposed  by  Oldham  et  al.  [13]  was 
followed  to  incorporate  the  monopodial  lung  structure  for  the  tracheobronchial  region.  Instead  of 
doubling  in  number  which  is  common  with  symmetric  dichotomous  lung  geometries,  the  number 
of  airways  per  generation  increased  by  1.5"  for  a  linear  increase  between  generation  5  (trachea 
being  considered  as  generation  0)  and  the  terminal  bronchioles  (1505  and  1584  for  B6C3F1  and 
BALB/c  mice,  respectively).  Airway  lengths  and  diameters  were  also  rescaled  to  lung  total 
capacity  for  consistency  with  other  lung  geometry  databases  in  MPPD.  Reported  airway 
measurements  of  Oldham  et  al.  [13]  and  Phalen  [14]  were  used  for  the  B6C3F1  lung  geometry. 
The  lung  geometry  of  BALB/c  mouse  was  from  the  measurements  of  Oldham  and  Phalen  [13] 
and  Oldham  and  Robinson  [15]. 

For  the  pulmonary  region,  Oldham  and  Robinson  [15]  for  the  BALB/c  mouse  and  Oldham  et  al. 
[13]  and  Phalen  [14]  for  the  B6C3F1  mouse  proposed  an  8-generation,  typical-path  dichotomous 
branching  alveolar  geometry.  The  alveolar  branching  structure  was  attached  to  the  distal  end  of 
the  terminal  bronchioles  in  the  tracheobronchial  tree  to  complete  the  lung  structure  for  BALB/c 
andB6C3Fl  mice. 

It  is  worth  noting  that  airway  parameters  had  to  be  rescaled  twice  prior  to  deposition  calculations 
in  each  strain  of  mouse:  first  by  (FRC/TLC)'^^ ,  where  FRC  is  the  functional  residual  capacity 
or  lung  volume  at  rest  and  TLC  is  the  total  lung  capacity,  to  adjust  airway  dimensions  to  rest 
conditions,  and  second  by  (1  -i-  Vj  /  2)^^^  ,  where  is  the  tidal  volume,  to  account  for  airway 
size  change  during  breathing.  The  rescaled  lung  geometry  was  used  in  the  calculations  of  particle 
deposition  in  the  lung. 

2.4  LUNG  VENTILATION  AND  PHYSIOLOGY  PARAMETERS 

Lung  ventilation  is  driven  by  the  difference  in  pressure  between  the  pleural  space  and  the  outside 
environment.  The  pressure  drop  due  to  the  lung  compliance  and  airway  resistance  will  produce 
airflow  in  lung  airways  during  inhalation  and  exhalation.  In  principle,  airflow  through  an  airway 
is  determined  from  the  solution  of  the  mass  and  momentum  balance  equations  for  the  entire  lung 
geometry.  Airflow  distribution  in  different  lobes  and  airways  of  the  lung  is  directly  related  to 
lung  expansion  and  contraction  rates  and  influenced  by  the  gravity-induced  intra-pleural  pressure 
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variation.  However,  lung  expansion  and  contraction  is  uniform  in  rodents  because  rodents  are 
typically  positioned  horizontally.  The  airflow  rate  entering  each  daughter  airway  of  an  airway 
bifurcation  under  such  circumstances  is  shown  for  uniform  lung  expansion  and  contraction  to  be 
proportional  to  the  distal  volume  to  each  daughter  airway  branch  [16].  Hence,  airflow  rate 
through  all  lung  airways  is  calculated  by  traversing  down  the  lung  tree  and  calculating  flow  rates 
in  daughter  branches  of  an  airway  bifurcation  from  the  parent  airflow  rate  and  distal-volume 
proportionality  of  daughter  branches. 

Various  physiological  parameters  such  as  ,  TLC ,  FRC ,  and  URT  volume  were  needed  to 
calculate  inhalation  and  exhalation  flow  rates  at  the  trachea  from  which  flow  rates  in  subsequent 
airways  could  be  determined.  The  measured  physiological  parameters  from  the  double  chamber 
plethysmograph  study  of  DeLorme  and  Moss  [17]  were  used  to  determine  default  breathing 
frequency,  inhalation  and  exhalation  breathing  fractions,  and  Vj  for  BALB/c  and  B6C3F1  mice. 
The  TLC  was  calculated  by  summing  up  the  individual  airway  volumes  and  alveolar  volumes  for 
the  BALB/c  and  B6C3F1  lung  geometries.  FRC  and  URT  volume  were  determined  by 
extrapolating  from  proportional  ratios  in  rats  ( FRC  /  TLC 30%  and  URT  /  TLC  '-'3%). 

2.5  MECHANISTIC  MODELOF  PARTICLE  DEPOSITION  IN  THE  LUNG 

The  particle  mass  conservation  or  convective-diffusion  equation  was  developed  by  accounting 
for  the  number  of  particles  that  entered,  exited,  deposited  and  remained  suspended  in  an  airway 
per  unit  time  [18-20].  It  was  assumed  that  particle  concentration  was  uniform  across  the  airway 
cross  section  but  varied  with  time  and  location  within  the  airway.  First,  particle  concentration 
reaching  the  trachea  was  calculated  by  including  the  filtering  effects  of  the  URT  by  Brownian 
diffusion  and  inertial  impaction.  Next,  the  convective-diffusion  equation  was  solved  in  the  lung 
geometry  of  B6C3F1  or  BALB/c  mouse  one  airway  at  a  time  starting  from  the  trachea  and 
traversing  down  the  lung  tree  to  find  particle  concentration  throughout  the  lung  during 
inhalation,  pause,  and  exhalation.  Particles  were  carried  through  lung  airways  by  the  airflow 
generated  by  a  uniform  expansion  and  contraction  of  lung  lobes  as  described  in  the  previous 
section.  Finally,  particle  deposition  in  each  airway  of  the  lung  was  calculated  by  integrating 
particle  flux  over  time  and  airway  volume.  Deposition  fraction  of  inhaled  particles  per  airway, 
lobe,  and  region  of  the  lung  was  found  by  dividing  the  calculated  deposition  by  the  mass  of 
inhaled  particles.  Particle  deposition  fraction  is  a  unique  property  for  a  given  lung  geometry  and 
breathing  parameters.  Local  or  regional  mass  deposited  or  number  of  particles  deposited  was 
calculated  by  multiplying  the  deposition  fraction  by  the  mass  or  number  of  particles  inhaled, 
respectively. 

2.6  SOFTWARE  IMPLEMENTATION 

We  modified  the  existing  MPPD  software  application  to  incorporate  the  B6C3F1  and  BALB/c 
mouse  respiratory  deposition  models  described  above.  First,  the  particle  inhalability,  nasal 
diffusion  efficiency,  and  nasal  impaction  efficiency  equations  were  implemented  in  the  MPPD 
code.  Next,  we  created  geometry  data  files  used  to  represent  the  trachebronchial  and  pulmonary 
airway  structure  in  both  strains  of  mouse.  The  files  contained  airway  length,  diameter,  branching 
and  gravity  angle,  and  alveolar  volume  (for  pulmonary  airways)  information  for  each  airway 
generation.  The  geometry  data  files  were  used  as  input  by  MPPD  to  form  the  tree  branching 
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structure  within  which  particle  transport  and  deposition  were  calculated.  Slight  modifications  to 
the  code  were  necessary  to  accommodate  the  monopodial  (instead  of  the  typical  dichotomous 
branching  structure)  nature  of  the  mouse  respiratory  tract.  Default  respiratory  tract  volume  and 
breathing  parameters  were  determined  for  the  mouse  model  based  on  available  measurements 
[17],  [21].  Particle  transport  and  deposition  within  the  mouse  lung  airways  was  calculated  in  the 
same  manner  as  other  animal  models  and  implemented  as  described  in  the  preceding  section. 

The  MPPD  graphical  user  interface  (GUI)  was  modified  to  accommodate  the  addition  of  the  two 
new  mouse  models.  The  B6C3F1  and  BALB/c  strains  are  now  among  the  options  in  the  input 
dialog  for  animal  geometry.  Once  either  mouse  model  is  chosen,  realistic  mouse  respiratory  tract 
volumes  (URT,  FRC,  TLC)  and  breathing  parameters  (Vx,  breathing  frequency,  and 
inspiratory/expiratory  fraction)  are  selected  as  default  starting  values  for  the  deposition  model. 
These  values  may  be  changed  at  the  user’s  discretion.  Once  these  values  are  entered  and  the 
inhaled  aerosol  characteristics  are  defined  by  the  user,  the  model  may  be  executed.  Deposition 
predictions  are  generated  by  region  (URT,  TB,  pulmonary,  LRT)  and  also  by  generation  number. 
Text  reports  containing  a  summary  of  the  model  predictions  and  tabulated  deposition  totals  as 
well  as  graphical  plots  of  deposition  predictions  may  be  generated  by  the  user.  These  time- 
stamped  reports  and  plots  also  contain  logs  of  the  input  values  used  for  the  predictions  and  may 
be  saved  and/or  printed  for  future  reference. 
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3.0  MODEL  PREDICTIONS 


Regional  deposition  in  the  LRT  was  ealeulated  for  B6C3F1  and  BALB/e  strains  of  miee  for 
whieh  lung  geometry  and  physiologieal  parameters  are  deseribed  above.  To  study  the  influenee 
of  lung  geometry  alone  on  partiele  deposition,  the  same  physiologieal  parameters  were  used  in 
the  eomputations,  whieh  ineluded  a  URT  volume  of  0.0315  cm^  [21],  lung  FRC  of  0.3  cm^  [17], 

Vt  of  0.2  cm^  [17],  average  breathing  frequeney  of  300  breaths/min,  and  0.4  for  inspiratory 
fraction  with  no  pause  between  inhalation  and  exhalation  [17]. 

Particle  deposition  in  the  mouse  lungs  is  suspected  to  vary  by  strain.  To  study  the  influence  of 
lung  geometry  on  particle  deposition,  particle  deposition  fractions  in  the  lungs  of  B6C3F1  and 
BALB/c  mice  were  calculated  via  endotracheal  breathing  (i.e.,  no  head  losses)  for  the  same  lung 
and  breathing  parameters  stated  above.  Predicted  deposition  fractions  are  given  in  Figure  3-1  for 
the  tracheobronchial  (TB)  and  pulmonary  (PUL)  regions  and  the  sum  (LRT)  of  deposition 
fractions  in  the  two  regions.  The  shapes  of  the  regional  deposition  fraction  curves  were  similar  in 
the  two  strains.  Particle  deposition  was  greater  in  the  BALB/c  mouse  lung  than  that  in  B6C3F1 
mouse.  TB  deposition  was  only  slightly  different  and  nearly  identical  for  0.1  pm  and  greater  size 
particles.  For  ultrafine  particles,  the  difference  in  TB  deposition  increased  with  decreasing 
particle  size  but  did  not  exceed  4%.  There  was  a  greater  inter-strain  difference  in  deposition 
fraction  in  the  pulmonary  region.  The  largest  difference  was  found  to  be  about  8%  for  0.03-  pm 
and  10-  pm  particle  sizes  in  the  pulmonary  region.  Overall,  there  was  a  maximum  9%  deposition 
difference  in  the  LRT  region,  which  occurred  for  0.02  pm  size  particles.  Thus,  adjustments  are 
necessary  to  correct  for  strain-specific  particle  deposition  in  the  mice  lungs. 


Figure  3-1.  Particle  deposition  in  the  lungs  of  mice  via  endotracheal  breathing. 
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The  calculated  deposition  fractions  of  particles  in  the  lung  are  plotted  against  generation  number 
in  Figure  3-2  for  the  case  of  nasal  breathing  to  examine  deposition  distribution  throughout  the 
LRT.  Deposition  fraction  decreased  for  both  strains  of  mice  with  increasing  particle  size  partly 
due  to  aerodynamic  properties  and  partly  due  to  the  reduction  in  respirability  with  increasing 
particle  size.  Deposition  of  1-  pm  particles  was  negligible  in  the  TB  region  and  peaked  in  distal 
generations  of  the  LRT.  Deposition  fractions  of  0.01-  and  0. 1-  pm  particles  increased  with  lung 
depth  (generation  number)  until  airborne  particles  were  depleted  from  the  inhaled  air  after  which 
deposition  fraction  dropped  to  zero.  The  results  indicated  that  to  reach  the  deepest  regions  of  the 
lung  with  significant  deposition,  particle  size  should  have  been  around  0.1  pm  .  Smaller  size 
particles  had  a  higher  upper  airway  deposition.  Larger  particles  demonstrate  reduced  respirability 
which  prevented  them  from  reaching  the  deep  lung.  It  is  interesting  to  note  that  no  significant 
impaction  deposition  occurred  in  the  first  few  upper  airways  of  the  lung  because  particle  sizes 
were  smaller  than  1  pm  .  The  results  also  confirmed  that  for  the  same  physiological  parameters, 
particle  deposition  was  higher  in  BALB/c  mice  than  in  B6C3F1  mice.  The  difference  was  most 
notable  in  the  alveolar  region,  which  has  the  greatest  deposition  fractions. 


Figure  3-2.  Deposition  fraction  of  various  size  particles  at  different  lung  depths. 

Regional  deposition  fractions  of  particle  with  sizes  between  0.01  and  10  pm  in  the  respiratory 
tracts  of  both  mouse  strains  are  given  in  Figure  3-3  via  nasal  breathing  for  the  same 
physiological  parameters  stated  above.  Comparison  of  predicted  deposition  fractions  in  the  LRT 
region  via  nasal  breathing  with  those  for  the  case  of  endotracheal  breathing  (Figure  3-1)  showed 
similar  patterns  for  submicrometer-sized  particles  but  drastic  differences  for  fine  and  coarse 
particles.  The  differences  were  caused  by  the  effects  of  particle  inhalability  and  nasal  deposition 
which  severely  reduced  the  number  of  particles  that  reached  the  LRT.  There  was  a  slight 
reduction  in  inhalability  and  some  deposition  in  the  nasal  passages  for  ultrafine  particles.  As  a 
result,  the  regional  deposition  curves  via  nasal  breathing  decreased  modestly  compared  with 
those  for  endotracheal  breathing.  Impaction  losses  increased  substantially  for  fine  and  coarse 
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particles  to  reach  100%  for  2  pm  particles  beyond  which  there  was  no  penetration  of  particles 
into  the  LRT.  Consequently,  the  regional  deposition  curves  for  nasal  breathing  dropped 
significantly. 


Figure  3-3.  Particle  deposition  in  the  lungs  of  B6C3F1  and  BALB/c  mice  via  nasal  breathing. 
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4.0  CONCLUSIONS 


A  particle  deposition  model  was  developed  for  the  lungs  of  BALB/c  and  B6C3F1  miee  and  used 
to  prediet  losses  of  particles  inhalable  to  humans.  Model  predietions  showed  that  only 
submierometer  particles  ean  penetrate  and  deposit  in  the  lungs  of  mice.  Ultra  fine  particle 
deposition  fraction  was  more  significant  in  the  pulmonary  region  than  in  the  TB  region.  There 
was  a  significant  deposition  of  large  particles  in  the  nasal  passages  of  miee.  The  deposition 
fraction  of  1  pm  or  larger  particles  was  only  about  a  few  percent.  Particles  larger  than  2  pm  were 
not  capable  of  penetrating  past  the  mouse  nasal  passages.  Comparison  of  model  predictions  for 
the  two  strains  of  mice  showed  that  while  the  deposition  patterns  were  similar,  there  were 
differences  in  the  predieted  deposition  fraction.  Thus,  strain-specific  deposition  models  should 
be  developed  for  mice  to  study  the  dose  and  biological  end  points. 
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6.0  DEFINITIONS,  ACRONYMS,  AND  ABBREVIATIONS 


ARA 

Applied  Research  Associates,  Inc 

B6C3F1 

Designation  for  a  specific  laboratory  mouse  strain 

BALB/c 

Designation  for  a  specific  laboratory  mouse  strain 

CDl 

Designation  for  a  specific  laboratory  mouse  strain 

DTRA 

Defense  Threat  Reduction  Agency 

FRC 

Functional  Residual  Capacity 

LRT 

Lower  Respiratory  Tract 

|Lim 

Micrometer  or  micron,  see  Conversion  Table 

MPPD 

Multiple-Path  Particle  Dosimetry  Model 

PUL 

Pulmonary 

Q 

Flow  rate 

TB 

T  racheobronchial 

TLC 

Total  Lung  Capacity 

URT 

Upper  Respiratory  Tract 
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